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The embryonic vertebrate heart tube develops an atrioventricular canal that divides the atrial 
and ventricular chambers, forms atrioventricular conduction tissue and organizes valve 
development. Here we assess the transcriptional mechanism underlying this localized 
differentiation process. We show that atrioventricular canal-specific enhancers are GATA- 
binding site-dependent and act as switches that repress gene activity in the chambers. We 
find that atrioventricular canal-specific gene loci are enriched in H3K27ac, a marker of active 
enhancers, in atrioventricular canal tissue and depleted in H3K27ac in chamber tissue. In the 
atrioventricular canal, Gata4 activates the enhancers in synergy with Bmp2/Smad signalling, 
leading to H3K27 acetylation. In contrast, in chambers, Gata4 cooperates with pan-cardiac 
Hdaci and Hdac2 and chamber-specific Heyl and Hey2, leading to H3K27 deacetylation 
and repression. We conclude that atrioventricular canal-specific enhancers are platforms 
integrating cardiac transcription factors, broadly active histone modification enzymes and 
localized co-factors to drive atrioventricular canal-specific gene activity. 
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Congenital heart defects affect ~ 1% of live births, and are 
found in up to 10% of spontaneously aborted fetuses^ 
Better insight into the transcriptional networks governing 
heart development would help to shed light on the origin and 
aetiology of congenital heart defects. A highly conserved step in 
the evolution of vertebrates is the division of the heart in atrial 
and ventricular chambers by the atrioventricular (AV) canal^. 
This step is crucial for heart structure and pump function. Defects 
in AV canal development account for a large fraction of 
congenital heart defects^"^. 

The initially formed heart tube consists of embryonic 
myocardium that proliferates slowly and conducts the electrical 
impulse at a slow rate^. While the tube elongates, specific regions 
differentiate and start to proliferate to form the expanding atria 
and ventricles. In contrast, the region in between the atrial and 
ventricular chambers does not expand but forms the AV canal. 
The cells of the AV canal preserve characteristics of the 
embryonic myocardium. The AV canal is required for septation 
and valve formation^, and forms pacemaker tissues that slowly 
conduct the impulse from atria to ventricles ensuring their 
sequential contraction^. During development, embryonic 
myocardial cells are thus confronted with a critical decision: 
differentiation into pacemaker-like cardiomyocytes of the AV 
canal or into working cardiomyocytes of the atrial and ventricular 
chambers. This decision is, at least in part, imposed by Bmp2- 
mediated signalling and the T-box transcription factors Tbx2 and 
Tbx3, which are selectively active in the AV canal where they 
suppress a default chamber gene programme^"^^. However, how 
AV canal-specific gene expression is regulated is not understood. 

An enhancer upstream of chicken Gata6 directs transgene 
expression throughout the primary heart tube and at later stages 
in the AV canal Furthermore, enhancers of Thx2 (ref. 13) 
and Gjd3 (Cx30.2)^'* and the promoter of Tnni3 (cTnl)^^ confer 
AV canal- restricted gene activity in transgenic mouse embryos. 
Thx2 was found to be activated by Bmp- signalling mediating 
Smads^^, whereas the activity of the Cx30.2 enhancer depends on 
the broadly expressed factors Tbx5 and Gata4 (ref. 14). 
Nevertheless, the mechanism for the AV canal specificity has 
remained unclear. Here, we investigated the regulatory 
mechanism underlying the AV canal- specific activity of these 
enhancers. We describe a mechanism by which GATA-dependent 
AV canal- specific enhancers mediate tissue- dependent gene 
activation and repression. GATA-binding elements in such 
enhancers recruit a Gata4/Smad4/HAT transcriptional 
activation complex in the AV canal and a Gata4/Heyl, 21 
HDAC transcriptional repression complex in the chambers. The 
former induces H3K27 acetylation that is associated with AV 
canal gene activation, the latter induces H3K27 deacetylation 
associated with repression of AV canal genes in the chamber 
myocardium. 

Results 

GATA-dependent AV canal enhancers act as chamber 
repressors. We first tested whether a mechanism of repression 
could be involved in the restriction of the activity of AV 
canal-specific enhancers to the AV canal. The promoter of the 
rat cardiac troponin T (cTnT, Tnnt2) gene drives efficient pan- 
cardiac expression in transgenic mice^^ (Fig. la,b). Previous 
analysis revealed that a tandem array of four copies of a 102 bp 
( — 1420/ — 1319) cGata6 enhancer module recapitulates AV 
canal specificity^^. This tandem array was coupled to the cTnT 
promoter and the activity was assayed in El 1.5 transgenic 
embryos. Strikingly, the cGata6 module restricts activity of the 
cTnT-lacZ transgene construct to the AV canal, indicating that it 
acts as repressor of the cTnT promoter in the chambers 



(Fig. lc,d). The cGata6 module does not appear to be 
conserved in mammals. Therefore, we next tested the conserved 
AV canal- specific mouse Tbx2 enhancer An array of four 
copies of the 380 bp enhancer likewise restricted cTn T-driven 
expression to the AV canal (Fig. le,f), indicating this was not a 
specific property of the cGata6 module, but a general property of 
AV canal- specific enhancers. We conclude that these AV canal 
enhancers act as tissue-specific switches, activating gene 
expression in the AV canal and repressing gene expression in 
the chambers. 

The AV canal- specific enhancers of cGata6, Cx30.2, Tbx2 and 
promoter of cTnl contain GATA-binding sites^^"^"*'^^ 
(Supplementary Fig. 1). Gata4 chromatin immunoprecipitation 
(Chip) -sequence data from adult hearts revealed that Gata4 
binds the AV canal- specific regulatory elements in vivo (Fig. li). 
This was confirmed by chromatin immunoprecipitation 
quantitative polymerase chain reaction using extracts prepared 
form AV canal or from the chambers (Supplementary Figs 2,3b). 
To assess whether the GATA sites are required for AV canal 
activity, we made use of a mouse heart organ culture and 
transfection assay to measure enhancer activities. We observed 
that mutation of the GATA-binding sites affected the activation 
of both the cGata6 and Tbx2 enhancer in AV canal explants 
(Supplementary Fig. 3d). 

Next, the GATA-binding sites were mutated in the cTnT 
promoter-coupled Tbx2 enhancer. We observed loss of repression 
of the cTnT promoter in the chambers in vivo (Fig. lg,h). These 
data revealed that the GATA sites are required for the repressive 
activity of this enhancer. To identif)^ additional binding elements 
involved in the activity of AV canal enhancers, we contrasted 
Gata4 ChIP sequencing data^^ and a set of AV canal-enriched 
genes^^. This analysis indicated that Gata4-binding sites close to 
or in AV canal- enriched genes are often associated with the 
presence of the binding site motif for the Sp family of 
transcription factors (Supplementary Fig. 4a). The 47 bp core 
module of the cGata6 enhancer^^ that contains a Sp-binding 
motif was mutated to only preserve the GATA-binding sites 
(cGata6muSP(GATA)2). Eight copies of the mutated module 
(cGata6muSP(GATA)i6) were coupled to the cTnT-promoter 
lacZ reporter. We found that the ability of the module with 
the mutated Sp site to repress transcription in the chamber 
myocardium was preserved, but with lower efficiency 
(Supplementary Fig. 4c). We conclude that the GATA sites are 
sufficient to mediate transcriptional repression in chamber 
myocardium. 



Assessment of H3K27 acetylation in AV canal and chambers. 

Gata4 itself is not localized, but is present in the nuclei of both 
AV canal and chamber myocardium in identical subnuclear 
localization patterns (Supplementary Fig. 3c). Therefore, we 
speculated that Gata4 cooperates with a localized factor that 
mediates gene activation in the AV canal or gene repression in 
the chambers. Previous work identified Gata4 as a regulator of 
histone modifications by physically interacting with Hats and 
Hdacs^^. Therefore, we tested whether AV canal- specific 
enhancers are differentially marked by acetylation of histone 
H3 lysine 27, a modification associated with active enhancers and 
with Hat p300 activity^ To this end, the genome-wide 
distribution of H3K27ac was assessed by ChlP-seq of purified 
embryonic AV canal cells and of chamber cells. We generated 
El 0.5 double transgenic embryos in which the expression of the 
green fluorescent protein (EGFP) (Tbx3BAC-GFP^) is restricted 
to the AV canal and of the red fluorescent protein Katushka 
(NppbBAC-Katushka^^) is restricted to the chamber myocardium 
(Fig. 2a,b). EGFP and Katushka-positive cells were purified by 
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Figure 1 | GATA sites are required to mediate repression of AV canal genes in chamber myocardium. Localization of transgene expression in E11.5 
mouse hearts. (a,b) The cTnTtransgene is expressed both in the AV canal myocardiunn and the channber myocardium. (c,d) cGato6-c7n7 transgenics show 
predominant expression in the AV canal, similar to the pattern of the 7bx2-c7n7 transgenes (e,f). (g,h) Mutation of the two GATA sites in the Tbx2 
regulatory region removes the repression in the chambers. Ratios denote proportion of embryos showing chamber repression. Scale bar represents 100 [xm. 
(la) left atrium; (ra) right atrium; (avc) AV canal; (Iv) left ventricle; (rv) right ventricle; (oft) outflow tract, (i) Browser view of Tbx2 gene locus with 
ChlP-seq profiles of Gata4 (ref. 18) (black) and p300 (ref. 24) (black) in heart, Smad4 (grey) in human ESC^^ and Hey2 (blue) in mouse ESC. The 380 bp 
Tbx2 (- 2700/ - 2300) enhancer region is depicted as a black box. 



fluorescence- activated cell sorting (Fig. 2c, Supplementary Fig. 5) 
and used for ChlP-seq (Fig. 2e). We asked whether DNA regions 
differentially enriched between the AV canal and chambers for 
H3K27ac were selectively associated with genes specifically 
expressed in the AV canal that were previously identified^^ 
(Supplementary Tables 1,2). The loci of the 50 genes whose 
expression is most significantly enriched in embryonic AV canal 
were also significantly enriched in AV canal- specific H3K27ac 
marks, indicating that activity in the AV canal is associated with 
H3K27 acetylation. The enrichment was validated in a small set of 
genes whose specific expression in AV canal myocardium has 
been established (Cacnalg, Cacna2d2, Kcna6, Msx2, Tbx2, Thx3, 
Idly Id2, Gjd3, Robol and Bmp2, Supplementary Tables 1,2). 
Furthermore, H3K27ac marks in loci harbouring genes expressed 
in both AV canal and chambers (including Gata4, Gata6, Nkx2-5, 



Tbx5, etc.) were distributed equally between AV canal and 
chambers. Chamber-specific H3K27ac marks were enriched in a 
set of chamber- specific gene loci {Hey2, Nppa, Nppb, Fig. 2e 
and ChlP-seq data set). This was confirmed by chromatin 
immunoprecipitation quantitative polymerase chain reaction 
from purified AV canal and chamber cells (Fig. 2f). For broadly 
expressed genes such as Gata4, our ChlP-seq data correlated with 
the H3K27ac E14.5 heart data of ENCODE^^ (Fig. 2e), whereas 
AV canal-specific H3K27ac signals were nearly undetectable in 
the ENCODE set as expected (Fig. 2e; E14.5 whole heart is 
composed of a very small fraction AV canal). 

The occupancy profile of Gata4 colocalized with regional 
H3K27 acetylation on AV canal enhancers (Fig. li. 
Supplementary Fig. 2). Gene ontology analysis of the gene sets 
co-occupied by Gata4 and either AV canal-specific H3K27ac or 
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Figure 2 | Repression of GATA-dependent AV canal enhancers in chamber myocardium is mediated by Hdacs and deacetyiation of histone H3K27. 

Whole-mount fluorescence nnicroscopy of double transgenic E10.5 ennbryo (a) and E10.5 heart (b). (c) Fluorescence-activated cell sorting profile of double 
AV canal-rcFP and channbers-Katushka transgenic E10.5 hearts. Scale bars represent 100 |im. (d) Chromatin immunoprecipitation quantitative polymerase 
chain reaction (ChlP-qPCR) assay with H3K27ac antibody using EGFP and Kat sorted cells. Bars represent mean±s.d. (n = 3). (e) Alignment of 
H3K27a ChlP-seq profiles from EGFP (green), Kat (red) sorted cells and E14.5 whole hearts^^ (black), (f) ChlP-qPCR assay with H3K27ac antibody using 
dissected AV canal/outflow tract and chamber myocardium extracts. Bars represent mean±s.e.m. (n = 3). Statistical test was conducted using the 
Student's two-tailed t-test (*P<0.05). (g) E10.5 cGata6-hsp68 transgenic hearts were cultured in the presence or absence of 50 nmol I ~^ TSA for 48 h. 
TSA-treated littermate hearts displayed pan-myocardial /acZ expression. In control medium, the reporter gene remained confined to the AV canal. Scale bar 
represents 100 (im. (h) E9.5 chambers-Katushka/AVcanal-EGFP transgenic hearts were cultured in the presence or absence of 50 nmol I ~^ TSA for 48 h, and 
expression of AV (Tbx2, Tbx3, EGFP) and chamber markers (Nppa, Nppb, Katushka) was analysed with real time RT-PCR. Data were normalized to HPRTand 
expressed as folds of increase over untreated samples, (i) Luciferase reporter assays on 102 bp cGataS, the 380 bp Tbx2, 660 bp Cx30.2 and the 356 bp cTnl 
genomic fragments. Cos-7 cells co-transfected with Gata4 expressing vector or not were treated in the absence or presence of 30 nmol l~^ TSA. Bars 
represent mean±s.e.m. (n = 6). Statistical test was conducted using the two- and three-way analysis of variance (*P<0.05 for Gata4 or TSA treatment 
versus control, *P<0.05 for Gata4 and TSA treatment versus Gata4 or TSA treatment). 
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chamber- Specific H3K27ac revealed that both sets were enriched 
for terms related to heart development (Supplementary Table 3). 
The Gata4-AV canal-specific H3K27ac set was enriched for terms 
related to cell metabolism. In contrast, the Gata4- chamber- 
specific H3K27ac set was enriched in genes linked to muscle 
differentiation and muscle cytoskeleton organization, which is in 
line with the notion that chamber myocytes acquire more mature 
cardiac muscle properties compared with the AV canal myocytes. 
These data indicate that in the two compartments different gene 
loci are both occupied by Gata4 and marked by H3K27ac. 

To assess whether the repression of AV canal enhancers in 
chamber myocardium is mediated by Hdac-mediated deacetyla- 
tion, explanted embryonic mouse hearts of the cGata6-lacZ line 
were treated with trichostatin-A (TSA) that blocks Hdac 
activity^^. When cultured in control medium, the reporter gene 
remained confined to the AV canal whereas TSA-treated 
littermate hearts displayed pan-myocardial lacZ expression 
(Fig. 2e). Expression of AV canal-specific Tbx2, Tbx3 and EGFP 
driven by Tbx3 regulatory regions increased in explanted E9.5 



embryonic hearts treated with TSA compared with 
dimethylsulphoxide- treated explants (Fig. 2h). We next assessed 
the activity of the 380 bp Tbx2 enhancer, the 660 bp Cx30.2 
enhancer and the 356 bp cTnl promoter in TSA-treated cells 
(Fig. 2i). Inhibition of Hdacs increased the lucif erase activity 
of the reporter constructs (Fig. 2i). These enhancers were 
synergistically and strongly induced by the combination of both 
TSA and Gata4. 



Smad4 and Gata4 cooperatively mediate H3K27ac. Bmp2 
signalling is selectively active in the AV canal myocardium^^. Its 
downstream effector Smad4 can physically interact with Gata4^^ 
and recruit the histone acetyltransferase p300 to specific 
chromatin loci^^"^^. We evaluated the transcriptional regulation 
of AV canal enhancers by Gata4 and Bmp2 signalling. The 
cGata6 enhancer (three tandem copies) showed a 7-fold 
activation by Gata4 in mouse C2C12 myoblast cells (Fig. 3a). 
The enhancer was not activated by Bmp2 alone, but was 



30 1 C2C12 



25 



^ 20 



^ 15 



2 10 



0 I rn PiPi f 
BMP2 300 ng mM 
Gata4 0.1 |ig 



600 
S 500 

1 400 
CO 300 

CD 
W 

2 200 
1 100 



Gata4 0.1 |ig 
Alk3 0.2 
Smadi 0.1 
Smacl4 0.1 



□ minP 

□ cGata6_minP 

□ Tbx2_minP 

■ Cx30.2_minP 

■ cTnLminP 



COS-7 

* 

I 



BMP2 E9.5 hearts 



_3.5 

§2.5 
o 

Z 2 
> 

^ 1.5 

CD 

< 1 

i 0.5 
^ 0 



www 



C\] CO Q_ 
X X li_ 

f ^ o 




Figure 3 | The Gata4-Smad-Hat complex induces H3K27 acetyiation and AV canal gene activation, (a) Luciferase reporter assays on the 102 bp cGata6 
enhancer. Plasmids encoding full-length Gata4 protein, the luciferase reporter construct with a nninimal promoter (nninP) or containing the cGata6 
enhancer with the same minimal promoter (cGata6_minP) were co-transfected in C2C12 cells. The cells were incubated with or without the recombinant 
human BMP2. Statistical test was conducted using two- and three-way analysis of variance (ANOVA) (mean ± s.e.m., n = 3, *P<0.05 for Gata4 or 
BMP2 treatment versus control, *P<0.05 for Gata4 and BMP2 treatment versus Gata4 or BI\/1P2 treatment), (b) Luciferase reporter assays on the 102 bp 
cGata6, the 380 bp Tbx2, 660 bp Cx30.2 and the 356 bp cTnl genomic fragments. Constructs were co-transfected with Gata4, Alk3CA, Smadi and Smad4 
expression vectors into Cos-7 cells. Luciferase activity was determined and normalized as fold over the reporter alone (mean ± s.e.m., n = 6, ^P<0.05 for 
Gata4 or BMP2 effectors versus control, *P<0.05 for Gata4 and BMP2 effectors versus Gata4 or BMP2 effectors, using two- and three-way ANOVA). 
(c) E9.5 chambers-Katushka/AV canal-EGFP transgenic hearts were cultured in the presence or absence of 200 ng ml ~^ Bmp2 for 48 h, and expression of 
AV canal (Tbx2, TbxS, EGFP) and chamber markers (Nppa, Nppb, Kat) was analysed with real time RT-PCR. Data were normalized to HPRTand expressed as 
folds of increase over untreated samples, (d) Untreated and Bmp2-treated E9.5 hearts were subjected to ChIP assay analyses with anti-H3K27ac antibody. 
Histograms represent enrichment for Tbx2, Cx30.2 and cTnl regulatory enhancers (mean ± s.e.m., n = 3, *P<0.05 using the Student's two-tailed Mest). 

(e) ChlP-qPCR analysis with Smad4 antibody on E10.5 embryo extracts from AV canal/outflow tract and chamber myocardium. Primers amplifying the BRE 
of the Id2 promoter serve as a positive ChIP control^^. Pulled-down DNA fragments were analysed with real-time PCR and enrichment indicates the 
ratio of ChlPed DNA to negative control regions, normalized for input DNA (mean ± s.e.m., n = 3, *P<0.05 using the Student's two-tailed t-test). 

(f) Cartoon depicting the overlapping expression patterns of Bmp2/Smad4 and Gata4 at stage E9.5-E10.5. 
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synergistically induced about 2 5 -fold by the combination of 
Bmp2 and Gata4. The 380 bp Tbx2 enhancer, the 660 bp Cx30.2 
enhancer and the 356 bp cTnl promoter were moderately induced 
by either Bmp effectors alone (constitutive active Alk3, Smadl 
and Smad4) or Gata4 alone. In contrast, co-transfection of both 
Bmp effectors and Gata4 produced a very strong synergistic 
induction in reporter activity (Fig. 3b). 

Next, in vitro cultured embryonic hearts of the chamber- 
Katushka/AV canal-EGFP transgenic mouse line (Fig. 2a,b) were 
treated with Bmp2 (Fig. 3c), which resulted in induction of 
endogenous AV canal genes {Tbx2, Tbx3) and reduction of 
chamber myocardium markers {Nppa, Nppb) (Fig. 3c). In vivo 
Chip assays demonstrated that the Bmp2- induced gene activity 
was associated with an enrichment of the activating H3K27ac 
epigenetic mark on AV canal-specific enhancers (Fig. 3d). 
Analysis of available Smad4 ChlP-seq data^^ from embryonic 
stem cells (ESC) revealed that Tbx2y Cx30.2 and cTnl regulatory 
regions are occupied by Smad4 in this context, suggesting that 
they can be co-occupied by Smad4 and Gata4 in the heart (Fig. li 
and Supplementary Fig. 2). Smad4 indeed occupies sites also 
occupied by H3K27ac in El 0.5 AV canal myocardium, but not in 
chamber myocardium (Fig. 3e). Consistently, although Smad4 is 
ubiquitously expressed at El 0.5, its localization is predominantly 
nuclear in the AV canaP^'^^. Taken together, the data indicate 
that AV canal-specific regulatory regions that harbour the GATA 
elements direct complex formation of Gata4, Smad4 and Hat, 
implicating these factors in AV canal- specific gene activity 
in vivo. 



AV canal enhancers are repressed by Hdacs and Heyl, -2. The 

AV canal enhancers were associated with H3K27ac in the AV 
canal, with non-acetylated H3K27 in chamber myocardium, and 
were activated by TSA treatment. Therefore, we tested whether 
the transcriptional repression was mediated by Hdacs. Nine 
Hdac family members are expressed in the developing heart^"*"^^. 
Using real time reverse transcriptase polymerase chain reaction 
(RT-PCR), we determined their expression levels in the 
EGFP+ AV canal cells and Katushka+ chamber cells 
fluorescence-activated cell sorting isolated from El 0.5 
Tbx3BAC-GFP;NppbBAC-Katushka mouse embryos. None of 
them was found to be significantly differentially expressed 
(Fig. 4b). We therefore speculated that a locally expressed 
Gata4 co-repressor might directly or indirectly recruit the broadly 
expressed Hdacs. Hdac2 can physically interact with Gata4 
(ref. 37). Furthermore, Hdacl and 2 are functionally redundant in 
the heart, and cardiac- specific inactivation of both Hdacl and 2 
resulted in neonatal lethality associated with cardiac 
arrhythmias^^. Therefore, we focused on Hdacl and 2. Both 
Hdacl and 2 could prevent the activation of AV canal enhancers 
by Gata4 and effectors of Bmp2 signalling in cell culture (Fig. 4a). 
Since dynamic changes in the subcellular localization of Hdacs 
can contribute to gene regulation^^, we evaluated the subcellular 
localization of Hdacl and 2 in El 0.5 hearts. Both Hdacl and 2 
were localized in the nucleus in both the AV canal and chamber 
myocardium (Fig. 4c,d). We therefore speculated that a locally 
expressed Gata4 co-repressor might directly or indirectly recruit 
the broadly expressed Hdacs. 

Several lines of evidence suggest that the Notch target genes 
Heyl and 2 could mediate the Gata4 dependent repression in the 
chambers. First, Heyl and 2 factors have been shown to act as 
transcriptional repressors by recruiting Hdacl (ref. 39). Second, 
although the regulatory sequences of cGata6, Cx30.2 and cTnl do 
not contain consensus Hey-binding sites (that is, E-box), it has 
been reported that Hey factors can mediate repression via the 
interaction with GATA factors^^'^^. Moreover, Heyl and 2 are 



selectively expressed in atrial and ventricular myocardium, 
respectively, but not in AV myocardium, and were shown to be 
involved in boundary formation between the AV canal and the 
chambers^^'^^. Analysis of Hey2 ChlP-seq data performed in 
mouse ESC revealed that the regulatory regions of Tbx2y Cx30.2 
and cTnl are occupied by Hey2 (Fig. li and Supplementary Fig. 2). 
When co-transfected, Heyl or Hey2 suppressed the Gata4-Smad- 
mediated activation of the AV canal enhancers in cell culture 
(Fig. 4e). To determine whether these AV canal enhancers were 
regulated by Notch target genes in vivo, cultured chambers from 
E9.5 Tbx3BAC-GFP; NppbBAC-Katushka mouse embryonic 
hearts were treated with the Notch inhibitor DAPT (Fig. 4f). 
The expression of marker of AV canal {Tbx2, Tbx3) increased 
(Fig. 4f). This was associated with an enrichment of the H3K27ac 
epigenetic mark on AV canal-specific enhancers (Fig. 4g). We 
conclude that these Gata4 responsive enhancers are suppressed in 
the chambers in response to Notch signalling (Fig. 4h). 

Compartment-specific co-factors function via GATA sites. To 

assess whether the GATA sites are necessary for recruiting acti- 
vators and repressors, we mutated the three GATA sites in the 
cGata6 enhancer and the two GATA sites in the Tbx2 enhancer 
(Fig. 5a,b). Mutation of the GATA sites affected the activation of 
the mutated enhancers by Gata4 (Fig. 5a). The increase in activity 
by addition of TSA was lost for the mutated enhancers (Fig. 5a). 
The co-activation of these mutated enhancers by Gata4 and Bmp2 
effectors was significantly reduced (Fig. 5b). To assess whether the 
GATA sites are sufficient for recruiting activators and repressors, 
eight copies of the mutated 47 bp cGata6 enhancer module (cGa- 
ta6muSP(GATA)i6) were coupled to a lucif erase reporter (Fig. 5c). 
Similarly to the wild- type 102 bp enhancer, the cGata6- 
muSP(GATA)i6 module was activated ~ 100-fold by Bmp2 
effectors and Gata4, indicating that the preserved GATA sites are 
sufficient for the synergistic activation by Gata4 and Bmp2 effec- 
tors. When Heyl or Hey2 where co-transfected, a significant 
inhibition of the stimulated reporter activity was observed. Simi- 
larly, HDACl and 2 were able to suppress the synergistic activation 
by Gata4 and Bmp effectors. These data are consistent with a 
model in which GATA sites serve as a platform to recruit co- 
factors that, in turn, mediate gene activation or repression, thus 
allowing GATA site -containing AV canal enhancers to act as tis- 
sue-specific transcriptional switches (Fig. 5d). 

Discussion 

We found that enhancers that mediate AV canal-specific gene 
activity are activated in the AV canal myocardium and repressed 
in chamber myocardium, thus acting as tissue-specific switches. 
Histone H3K27 acetylation marks active enhancers^^ Although 
the pattern of this modification has been assessed genome-wide in 
heart tissue by ChlP-sequencing^"*, the data do not provide cell 
type- or cardiac compartment- specific information (whole hearts 
were analysed, the heart contains many fibroblasts and other 
cell types). We assessed the genome- wide H3K27ac pattern in 
isolated embryonic AV canal and in chamber cardiomyocytes, 
respectively, and found prominent compartment- specific 
H3K27ac patterns. Specifically, we observed that most of the 
estabHshed AV canal-specific loci were strictly associated with 
H3K27ac in the EGFP-labelled AV canal cells. Conversely, except 
for firmly established chamber loci such as Nppa, Nppb and 
Hey2, we found very few loci that were strictly enriched in 
H3K27ac in the Katushka-labelled chamber cells. Therefore, our 
data implicate this epigenetic signature in the regulation of 
embryonic AV canal specificity. 

The AV canal enhancers appeared to function through GATA 
factor-binding sites that are required and largely sufficient for 
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Figure 4 | The Gata4/Hey/Hdac complex induces H3K27 deacetyiation and AV canal gene repression, (a) Luciferase reporter assays on 102 bp cGata6, 
the 380 bp Tbx2, 660 bp Cx30.2 and the 356 bp cTnl genomic fragments. Constructs were co-transfected with Gata4, Alk3CA, Smadl and Smad4, 
HDAC1 and 2 expression vectors into Cos-7 cells (mean ± s.e.m., n = 6, *P<0.05 for Gata4, BMP2 effectors and HDAC1, 2 vs Gata4 and BMP2 co- 
transfection using two- and three-way analysis of variance (ANOVA)). (b) qPCR analysis of mRNA isolated from AV canal-EGFP and chambers-Katushka 
sorted cells. Values are normalized against EGFP population. Immunohistochemistry for Hdaci (c) and 2 (d) proteins on sections of E10.5 hearts. Scale bars 
represent 100 |im. (e) Luciferase reporter assays on AV canal enhancers with expressing vectors encoding for Gata4, Bmp2 effectors and the notch 
target genes Heyl and 2 (mean ± s.e.m., n = 6, *P<0.05 for Gata4, BMP2 effectors and Heyl, 2 versus Gata4 and BMP2 effectors co-transfection, using 
two- and three-way ANOVA). (f) E9.5 chambers-Katushka/ AV canal-EGFP transgenic hearts were cultured in the presence or absence of 50|imol l~^ 
y-secretase inhibitor (DAPT) for 48 h, and expression of AV canal and chamber markers was analysed with real time RT-PCR. Data were normalized to 
HPRT and expressed as folds of increase over untreated samples, (g) Untreated and DAPT-treated E9.5 hearts were subjected to ChIP assay analyses 
with anti-H3K27ac antibody. Pulled-down DNA fragments were analysed with RT-PCR and enrichment indicates the ratio of ChlPed DNA to negative 
control regions, normalized for input DNA. (h) Cartoon depicting the overlapping expressions of Heyl, 2 and Gata4 at stage E9.5-E10.5. 



both activation and repression functions. Mutation of GATA- module in which only GATA sites were maintained was sufficient 
binding sites with confirmed Gata4 occupancy, abolished to suppress activity in the chambers. We therefore identified a 
activating and repressive activities, and a minimal enhancer class of AV canal- specific enhancers essentially composed of 
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Figure 5 | GATA sites serve as a scaffold to recruit functionally distinct complexes. (a,b) \n vitro reporter assays with the cGata6 (102 bp) and Tbx2 
(380 bp) luciferase constructs containing mutations for GATA sites (cGata6muGATA, Tbx2muGATA). Cos-7 cells co-transfected with Gata4, Alk3CA, 
Smadi and 4, Heyl, 2, HDACl and 2 expressing vectors or not were treated in the absence or presence of 30 nnnol I ~ ^ TSA (nnean ± s.e.m., n = 3, ^P<0.05 
for Gata4 or TSA/BMP2 treatment versus Gata, TSA or BMP2 alone, *P<0.05 for mutated versus wild-type enhancers using two- and three-way analysis 
of variance (ANOVA)). (c) In vitro reporter assays on 102 bp cGata6 and mutated 47 bp cGata6 enhancer (cGata6muSP). Cos-7 cells were co-transfected 
with Gata4 and Gata4 co-activators and co-repressors (mean ± s.e.m.; n = 3; *P<0.05 for HDACl, 2 or Heyl, 2 versus Gata4 and BMP2 effectors; 
*P<0.05 for HDACl, 2 and Heyl, 2 versus HDACl, 2. ^P<0.05 for HDACl, 2 and Heyl, 2 versus Heyl, 2 using three-way ANOVA. (d) GATA-binding 
site-enhancers recruit a Gata4/Smad4/Hat transcriptional activation complex in the AV canal and a Gata4/Hey1,2/Hdac transcriptional repression 
complex in the chambers, which coordinately establish AV canal gene specificity. 



GATA-binding sites and marked by AV canal-specific H3K27ac. 
Although the genome contains thousands of occupied GATA 
sites, it is not clear why only a fraction of these sites have AV 
canal-specific properties. We therefore assume that AV enhancers 
contain unidentified binding sites for additional factors that 
specif)^ GATA site- dependent complex recruitment to drive AV 
canal-specific gene expression. 

Gata4 and Gata6 are essential for cardiogenesis. They are 
co-expressed in the myocardium and function redundantly"*^. 
Focusing on Gata4, we observed its presence in the nuclei of both 
AV canal and chambers coinciding with enhancer occupancy in 
both these compartments. Therefore, we explored the possibility 
of the presence of compartment- specific co-factors to try and 
explain the spatial restriction of AV canal enhancer activity. 

8 



Bmp2, an AV canal-specific Bmp ligand, is required for AV 
canal-specific gene expression, although the mechanism of 
enhancer activation has not been established *^'^^''*'*''*^. We 
observed that Smads (Smad4, Smadl/5), downstream effectors 
of Bmp2 (ref. 46), strongly activate the enhancers in a Gata4 and 
GATA-binding site-dependent manner. The in vivo H3K27ac 
Chip data support a mechanism in which the activating Gata4/ 
Smad complex recruits Hats such as p300 to a specific subset of 
AV regulatory regions, thereby inducing H3K27 acetylation and 
activation of AV canal genes. 

We identified a novel repressive mode of AV canal enhancers 
in the chamber myocardium, which restricts their activity to the 
AV canal. This observation suggests that, in addition to the 
localized Bmp2-mediated activation mechanism, an additional 
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layer of regulation is involved in limiting enhancer activity to the 
AV canal. One of the AV canal enhancers studied in this 
respect regulates Tbx2, a T-box factor critical for AV 
canal development. This factor effectively blocks chamber 
differentiation^^''*^"'*^, indicating that delimiting its activity to 
the AV canal is critical for establishing correctly localized 
AV canal development. The transcription factors that activate 
cardiac and AV canal enhancers are typically expressed in both 
AV canal and chamber myocardium (Supplementary Table 1). 
Furthermore, non-Bmp2-mediated Bmp- signalling is present in 
the chamber myocardium. Loss of function of Bmp6, -7 and -10, 
BmpRlA (Alk3) and Smad4 results in chamber and septal 
defects^^"^"*. The repressive activity in chamber myocardium we 
identified may therefore be required to suppress activation of AV 
canal- specific enhancers in the chambers by these spatially less 
unrestricted activators. 

Our results indicate that in chamber myocardium, Gata4 
cooperates with Hdacl and 2 and induces H3K27 deacetylation of 
AV canal gene loci, providing the context for repressing gene 
activity. Whether, and to what extent, this repression is chamber- 
specific, is not clear. None of the Hdacs examined showed any 
specificity for chamber myocardium. Furthermore, Hdacl and 2 
were localized in the nuclei of both AV canal and chamber 
myocardium (Fig. 4), and therefore are available for their 
deacetylation and repression functions. Our data suggest that 
Heyl and Hey2, expressed in atrial chamber and ventricular 
chamber myocardium, respectively"*^ and able to interact with 
Gata4 (ref. 40) and thus provide chamber-specificity to the Gata- 
Hdac activity. Furthermore, we identified a common, AV canal 
cell type-specific chromatin signature for AV canal-specific loci, 
which includes H3K27ac and co-occupancy of Gata4, Smad4 and 
Hey2 proteins. Previous work demonstrated that deficiency of 
both Heyl and 2 results in moderate expansion of Tbx2 expression 
in the chambers"**'^^, suggesting that these factors contribute, but 
are not solely responsible for chamber repression by Gata4-Hdac. 

Previous work revealed that chromatin remodelling enzymes 
regulate Gata4 target genes by post-translational modifications of 
Gata4 (ref. 20). Gata4 was found to interact with Hdac2 and Hopx, 
a homeodomain-only protein required for heart development^^. In 
the adult heart, Hopx is enriched in and required for the function 
of the AV conduction system^^. Hdac2 -mediated deacetylation of 
Gata4 may result in repression of Gata4 target genes in the AV 
canal and the derivative tissues. The PRC2 complex, a critical 
repressor required for normal lineage- specific differentiation, 
methylates Gata4, attenuating its transcriptional activity by 
reducing its interaction with and acetylation by p300 (ref. 57). 
Currently, we do not have data regarding the compartment 
specificity of such post-translational modification activities that 
could contribute to the spatial activity of the GATA-binding site- 
dependent AV canal enhancers. However, we hypothesize that the 
Gata4-Smad-HAT complex present in the AV canal would 
effectively compete with this mechanism and result in the 
activation of GAT A site- dependent AV canal enhancers (Fig. 5d). 

The AV canal myocardium is a central component of the 
cardiac conduction system in vertebrates^^'^^. It has pacemaker 
properties and its gene expression profile is distinct from 
that of the chamber myocardium^^'^^. These properties are 
reminiscent of those of the embryonic heart tube, implicating 
that the mechanism dictating AV canal specificity is also active in 
the embryonic heart tube prior to the development of the 
chambers, but is repressed in the differentiating chamber 
cardiomyocytes^^'^^. Indeed, the cGata6 enhancer in mouse and 
hmp4 and thx2h in zebrafish are active in the embryonic heart 
tube and maintain their primary, AV canal specificity by 
preventing their activity in the developing chambers, strikingly 
illustrating this process '^^'^^ Confining this gene programme to 



the AV canal may be essential for the electrophysiological 
integrity of the heart. Mutations in Gata4 and Gata6 are 
responsible for congenital heart defects including AV canal 
defects and arrhythmias in humans and mouse models*^'^^'^^'^^, 
implicating their function in both chamber and conduction 
system tissues. Inactivation of both Hdacl and 2 in the heart 
resulted in neonatal lethality associated with cardiac 
arrhythmias^^. Interestingly, expression of conduction system 
genes Tnnil, Cacna2d2 and Cacnalh^^ was upregulated in the 
heart chambers of Hdacl and 2 double mutant mice^^. Moreover, 
a regulatory region of the pacemaker channel Hcn4, active in 
pacemaker cells of the conduction system, was ectopically 
activated in the chambers by Hdac inhibition^^. The broadly 
expressed factor Proxl recruits Hdac3 and negatively regulates 
Nkx2.5 to ensure the functional and structural integrity of the 
pacemaker fijnction in adult mice^^. The Proxl-Hdac3-Nkx2.5 
pathway appeared to be Hdac2 independent in cells^^. This 
highlights the important role of Hdacs and chromatin -mediated 
mechanisms in the patterning of the AV conduction system. Our 
data suggest that Gata4-Hdac-mediated deacetylation of H3K27 
serves as a mode of confinement of the conduction system gene 
programme, preventing its activity in the chambers. Dysfunction 
of this mechanism may result in misexpression of pacemaker 
genes in the chambers, resulting in arrhythmias. 

Methods 

Ethics statement. All animal work conformed to ethical guidelines and was 
approved by the animal ethical committee of the Academic Medical Center of 
Amsterdam. 

Transgenic reporter assay. Each multimer of AV canal enhancer fragments was 
built in a geometric progression by first cloning a unit length PGR fragment, or 
annealed oligonucleotides in the case of the 47 bp and 102 bp cGata6 fragments, 
with Sall/Xhol ends into a similarly digested plasmid. Tetramers of the 380 bp 
{Tbx2), 102 bp and 47 bp {cGata6) fragments were released with Sall/Xhol, blunt- 
ended with Klenow and cloned into the Smal site upstream of the promoter in the 
cTnT-lacZ construct. Tbx2muGATA-cTnT is identical to the Tbx2-cTnT construct, 
with the exception of 4-bp substitutional mutations of the GATA sites located at 
position — 2312 and position — 2345 of the Tbx2 enhancer region. (GATA, GAAG 
GATATTGT, GGTGGATAAAGA; GATAmu, GAAGTTTTTTGT, GGTGTTTT 
AAGA). The mutations were generated using the QuickGhangeTM Site-Directed 
Mutagenesis kit (Stratagene). Tandemly repeated copies of a 102 bp cGata6 
enhancer fragment were cloned into the hsp68-lacZ reporter vector as previously 
described. DNA was injected into the pronucleus of 0.5 -day- old fertilized FVB/N 
eggs, which were subsequently transferred into the oviducts of GD-1 pseudo- 
pregnant foster females (Gyagen). At least 200 injections were performed per 
construct. Embryos were harvested and stained with X-gal to detect lacZ activity, 
and yolk sacs were processed for PGR genotyping. Only patterns that were 
observed in at least three different embryos resulting from independent transgenic 
integration events of the same construct were considered reproducible. A trans- 
genic mouse line harbouring a BAG of the Nppb locus modified with the red 
fluorescent protein Katushka^^ was bred with a transgenic mouse line harbouring a 
BAG of the Tbx3 locus modified with the EGFP^^. 



Embryonic explant culture. Mouse E9.5 embryonic hearts explants, dissected in 
Dulbecco's modified Eagle's medium (DMEM), were cultured in DMEM supple- 
mented with 5% fetal calf serum and penicillin/streptomycin for 48 h at 37 °G in 
5% GO2. Trichostatin-A (Sigma), BMP2 (R&D Systems), Inhibitor IX (Galbio- 
chem) or dimethylsulphoxide alone was added. 

Cell sorting. El 0.5 double transgenic hearts were dissected, pooled and digested 
with 0.25% Trypsin/EDTA (Invitrogen), neutralized in DMEM (Invitrogen) 
containing 5% FBS and lOmmoll" ^ HEPES (Invitrogen), rinsed and resuspended 
in phosphate-buffered saline (PBS) and passed through a 70-mm nylon cell strainer 
(Falcon). Samples were sorted on a FacsAria flow cyto meter (BD) using Influx 
software. Samples were gated to exclude debris and cell clumps. The number of 
El 0.5 AV canal GFP cells and chambers- Katushka cells per embryo obtained were 
typically 600-900 and 2,400-3,600, respectively. Fluorescent cells were collected 
into PBS and processed for RNA extraction or GhlP-seq. 

RNA isolation and real time RT-PCR. Total RNA was isolated from treated 
explants and sorted cells with NucleoSpin RNA XS (Bioke) following the protocol 
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of the manufacturer. cDNA was generated using Superscript II (Invitrogen) and 
OligodT as primers. Expression level of different genes was assessed with quanti- 
tative RT-PCR using the LightCycler Real-Time PGR system (Roche Diagnostics) 
and the primers described in Supplementary Table 4. Values were normalized to 
HPRT expression levels. 

Chromatin immunoprecipitation. E10.5 embryos were dissected, and the primary 
myocardium (AV canal/outflow tract) and atria and ventricle regions were separately 
collected in cold PBS. ChIP was performed with a True MicroChIP kit (Diagenode) 
according to the manufacturer's protocol. Antibodies used for ChIP were anti-Gata4 
(2 |ig, Santa Cruz Biotechnology), anti-Smad4 (2 |ig, Santa Cruz Biotechnology) and 
anti-H3K27ac (0.5 |ig. Abeam ab4729). Chipped DNA were subjected to library 
preparation using the 5500 Series SOLID Systems preparation kit (Applied Biosys- 
tems) according to manufactures recommendations and sequenced using the 5500 
wildfire sequencer (SOLID). Quantitative PCR was performed on a Roche LightCycler 
480 System using Sybr Green detection. Fold enrichment indicates the ratio of ChlPed 
DNA to negative control regions, normalized for input DNA. Primer sequences are 
provided in Supplementary Information. 

Hey2 ChIP was performed in CM7/1 mouse ESC with doxycycline-inducible 
Flag-Hey2 expression^^. 10^8 cells were induced with 50ng|a.l~ doxycycline 
for 48 h and ChIP was performed with an anti-Flag-M2 antibody (10|ig 
Sigma- Aldrich). Chipped DNA were subjected to library preparation using the 
NEBNext ChlP-seq sample preparation kit (New England Biolabs) according to 
manufacturer's recommendations and sequenced on a GAIIx platform (Illumina). 

X-gal staining and immunohistochemistry. For X-gal staining, embryonic hearts 
were dissected in PBS and then fixed for lOmin with 4% (wt/vol) paraformalde- 
hyde in PBS. Embryos were washed twice with PBS and then stained with X-gal 
(Invitrogen) for 4-16 h at 37 °C. After developing the blue colour, embryos were 
rinsed in PBS and post-fixed in 4% (wt/vol) paraformaldehyde. Fluorescent 
immunohistochemistry on sections was carried out with anti-Gata4 (1:50, Santa 
Cruz Biotechnology), anti-cTnl (Millipore), anti-Hdacl (1:200, Biovision) and 
anti-Hdac2 (1:200, Biovision) antibodies. Images were obtained with Leica DFC320 
and Leica DM6000. 

In vitro and ex vivo reporter assays. Tandemly repeated AV canal enhancers 
were cloned in pGL4-2.6-luciferase construct that contains a minimal promoter 
(Promega). Luciferase reporter constructs were co-transfected with expression 
constructs for human Gata4 (pcDNA3^^), mouse Smadl (pcDNA3^^), Smad4 
(pcDNA3^^), constitutively active Bmp-receptor Alk3 (Alk3CA) 
(pCS2.BmpRla.CA), Heyl, 2, HDACl and 2. Full-length mouse FLAG- Heyl, 2 
were kindly provided by Dr Esther Creemers (Academic Medical Center, 
Amsterdam, the Netherlands). Full-length human HDAC2 was kindly provided by 
Dr John Epstein. Full-length human HDACl was purchased from Open Biosystem. 
Constructs were transfected into Cos-7 cells with the PEI transfection reagent. Cell 
extracts and luciferase assays were performed following the protocol of the 
manufacturer (Promega). The mean luciferase activities and standard deviations 
(s.d.) were plotted as fold activation when compared with the empty expression 
plasmid. For DNA transfection of AV canal explants, 50 ng of firefly luciferase 
reporter vector was mixed with 80 ng |ig~ ^ of CMV |3-galactosidase plasmid in 5 |il 
of OPTIMEM (Invitrogen); 0.5 |il of Lipofectamine 2000 (Invitrogen) was added to 
5 |il of OPTIMEM. After 5 min, the DNA and Lipofectamine solutions were mixed 
together. After incubation for 20 min, the transfection reagent mixture (10|il) was 
transferred into a well of a 3 84- well plate containing the explant in medium. 
Transfection assays were performed in quadruplo and the efficiency was assessed 
by the expression of a CMV-EGFP reporter. Explants were harvested 48 h after 
transfection and processed for chemiluminescence activity using the Dual-Light 
luminescent reporter gene assay (Invitrogen). 

Processing of ChlP-sequencing data and statistical analysis. ChIP- seq data for 
H3K27ac in chamber- Katushka/ A V canal- EGFP- sorted cells chamber- Katushka/ 
AV canal-EGFP was processed using an in-house peak calling algorithm dubbed 
OccuPeak (van Duijvenboden et al, submitted). Using an identical threshold 
(P = 0.01), peak calling resulted in 87,572 and 67,092 H3K27ac peaks in AV canal 
and chamber tissue, respectively. 

Microarray data^^ was used to determine a top-50 of AV canal- enriched genes 
(Supplementary Tables 1,2). In addition, a custom list of firmly established AV 
canal- restricted genes was used (Supplementary Tables 1,2). The coordinates of 
these genes were imported from the UCSC gene database. Gene regulation by 
H3K27 acetylation is expected to occur on both promoters and distal enhancers. 
Therefore, we defined gene loci as UCSC genes plus flanking regions of 50 Kb on 
either side. Peaks were considered to occupy a gene locus if they overlap it with 
1 bp or more. Overlap was established using the Galaxy intersect tool (https:// 
www.main.g2.bx.psu.edu/) . 

Compared with the chamber H3K27ac ChlP-seq data, the AV canal data 
contained 1.4 times the number of mapped reads. Therefore, we assumed that the 
difference in total number of ChlP-seq peaks called between AV canal- (87,572) 
and chamber tissue (67,092) does not reflect biology, but is caused by a difference 
in peak calling power between the data sets. To correct for this difference, we 



determined the number of peaks from both data sets as present in 5,541 random 
gene loci. This resulted in a 0 hypothesis in which AV canal peaks are 1.08 times 
more common than chamber peaks (Supplementary Table 2). We assumed a 
normal distribution of H3K27ac ChlP-seq peaks across gene loci. The 1^ test was 
used to analyse the frequency of H3K27ac peaks according to tissue type and gene 
cluster. Statistical significance was set at P<0.05. 

Occurrence of specific motifs among the Gata4 peaks was performed with 
MEME^^. Biological process terms were identified using the DAVID informatics 
resource gene ontology analysis software^^. 

Statistical analysis. For gene expression analysis and ChIP assay, P-values were 
generated using the Student's two-tailed Mest. For luciferase assays analysis of 
variance was used. Statistical significance was considered for P<0.05. 
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